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ABSTRACT 



Context. Models of chemical evolution of galaxies including the dust are nowadays required to decipher the high-z 

universe. In a series of three papers we have tackled the problem and set a modern chemical evolution model in 

which the interstellar medium (ISM) is made of gas and dust. In the first paper (Piovan et al., 2011a) we revised the 

condensation coefficients for the elements that typically are present in the dust. In the second paper (Piovan et al., 

2011b) we have implemented the dust into the Padova model with infall and radial flows and tested it against the 

observational data for the Solar Neighbourhood (SoNe). In this paper (the third of the series) we extend it to the whole 

Disk of the Milky Way (MW). 

Aims. The Disk is used as a laboratory to analyze the spatial and temporal behaviour of (i) several dust grain families 

with the aid of which we can describe the ISM dust in a simple way, (ii) the abundances in the gas, dust, and total ISM 

of the elements present in the dust emitted by stars and the dust grains grown by accretion in the ISM, and finally (iii) 

the depletion of the same elements. 

Methods. The temporal evolution of the dust and gas across the Galactic Disk is calculated under the effect of radial 

flows and a Bar in the central regions. The gradients of the abundances of C, N, O, Mg, Si, S, Ca and Fe in gas and 

dust across the Disk are derived as a function of time. 

Results. The theoretical gradients nicely reproduce those derived from Cepheids, OB stars, Red Giants and HII regions. 

This provides the backbone for the companion processes of dust formation and evolution across the Galactic Disk. We 

examine in detail the contributions to dust by AGB stars, SNae and grain accretion in the ISM at different galacto- 

centric distances across the Disk. Furthermore, we examine the variation of the ratio between silicates and carbonaceous 

grains with time and position in the Disk. Finally, some hints about the depletion of the elements in regions of high 

and low SFR (inner and outer Disk) are presented. 

Conclusions. The results obtained for the Disk of the MW make it possible to extend the model we have developed to 

other astrophysical situations such as galaxies of different morphological types, like ellipticals or star-busters, or even 

to a different theoretical modelling, like the chemo-dynamical N-Body simulations. 

Key words. Galaxies - Dust; Galaxies - Spirals; Galaxies - Milky Way 



1. Introduction 



New generation, powerful telescopes and space instru- 
mentation have unveiled a high-z universe surpris- 
ingly obscured by large amounts of dust, and con- 
sequently spurred a new generation of chemical and 
spectral galaxy models in which the formation, pres- 
ence and ev olution of dust i s take n into accoun t (Se e 
for instance 



Zhu kovska et all 1 2008: Calu ra et all [2 008: 



Valiante et al. 



2009; Pirjino et al. 2011; Gall et al. 2011a,b 



Mattssonll20lU I Valiante et al J 120111; lAcharvva et al.l 12011 



Kemper et al.l 1201 it iDwek fe Cherchnefj 120111 to mention 
the most recent ones). Given the interest in understanding 
the properties of the early universe, most of these models 
are tuned toward QSOs or star-burst galaxies. They try to 
get clues about the role of dust in the early stages and to 
estimate how much dust can be produced both by stars 
and other processes in the interstellar medium (ISM) on 
a given timescale in order to match the amo u nt of dust 
observed in high-z objects dGall et al.l I2011allb1 ; iMattssonl 
1201 It iDwekfe Cherchnefj [20ilD . However, because of the 



huge distances, the information on dust to our disposal for 
these extremely far objects comes from (i) the extinction 
curve, that suggests clues about dust composition by the 
comparison with the one for galaxies in the local such as 
the SMC, LMC, and M; (ii) estimates of the dust mass 
from the fluxes observed in the pass-bands to our disposal. 
Other pieces of information are typically related to the 
upper limit to the age of the target, its metallicity and star 
formation rate (SFR) and, finally, the total masses in form 
of stars, molecular g as (if CO obse r vation s are available), 
and Dark Matter (|Valiante et al.l 1201 1[ ). Since we are 
observing the integrated properties of these galaxies, the 
history of chemical enrichment for individual elements 
and different distances from the centre of the galaxy are 
hardly available. A detailed information about the way in 
which the various elements behave is however required if 
we want to cross compare data with theoretical modelling 
of dust formation/evolution, because it would correlate 
the processes affecting the dust with the gas composition 
existing in that part of the galaxy we intend to model 
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(jZhukovska et all 120081: iPiovan et al.ll2011bl ). This indeed 
would be the typical case in which the relative proportion 
of silic ates /carbonaceous grains and iron dust develop with 
time ()Dwekl 120051: IZhukovska fc Gaill l2009t ) in different 
regions of a galaxy, thus getting hints on the dust mixture. 
For its proximity and hence richness of data, the best 
possible laboratory is the SN. Another useful laboratory 
in which we can trace the evolution of the element abun- 
dances by means of the observations of different types of 
objects is the Disk of the MW. Nevertheless, till now, only 
an handful of models have considered the presence of dust 
in the Galactic Disk and its evoluti on in space and time 
(|Dweklll998l IZhukovska fc Gailll2009lh 

Basing on the observational data we have referred 
to above, in this study we plan to examine the forma- 
tion/evolution of dust and depletion of the elements along 
the MW Disk, trying to satisfy the constraints derived 
from the radial gradients in matter and abundances of 
various chemical elements as measured in different types 
of sources. This paper is the last one of a series of three 
in which we analyzed the formation and evolution of 
dust in the MW. It is worth briefly summarizing the 
aims and results of the first two papers, and also briefly 
comment on the gradual develop ment of the chemical 
model. In the first paper (|Piovan et al.l l2011af) we ex- 
amined and presented a compilation of condensation 
efficiencies to be adopted to st udy the dust conten t of 
the gas. In the second paper (jPiovan et al.l l2011btl we 
apply our condensation coefficients to a detailed chem- 
ical model of the Solar Neighborhood and reproduced 
the local depletions of the chemical elements together 
with other observational data getting clues on the effect 
induced by important parameters such as the IMF, the 
star formation law, and the accretion model adopted to 
simulate the grain growth in the ISM. In addition to 
this, the same analysis has been made for the inner and 
outer regions of the MW Disk, trying to understand how 
they would behave in a high/low star forming environment. 

The plan of the present paper is as follows. In Sect. 
[21 we summarize the main features o f the chemical mode l, 
more details can be also found in (jPiovan et al.l 1201 lbt> . 
In particular, we cast the equations governing the model 
in presence of radial flow of matter and dust. In Sect. [3] 
we present the equations for the radial flows in presence 
of dust, gas, and ISM (gas+dust). In Sect. 13.11 we discuss 
the boundary conditions and the final equation for the in- 
nermost shell (Sect. I3.1.1J) . the generic shell (Sect. 13.1. 21) . 
and the outermost shell (Sect. l3~l~3l) . Some computational 
details are discussed in Sect. 13.21 In Sect. 2] we present a 
simple algorithm based on the Artificial Neural Networks 
(ANNs) to predict the amount of molecular gas (cold re- 
gions where we assume that dust accretion takes place) in 
the different rings of the Disk taking into account the infor- 
mation about SFR and surface mass density. In Sect. [5] we 
compare the results of the chemical models, whose proper- 
ties for the Solar Vi cinity have been already discussed in 
I Piovan et al.l (|2011bt ). extended now the whole Disk. The 
theoretical predictions are compared with the observational 
data for the gradients in gas content, chemical abundances, 
and abundance ratios. Furthermore, we present and discuss 
the radial and time evolution of many dust-related quan- 
tities. In particular, we analyze the radial behaviour of (i) 



some typical grain families; (ii) the elements involved in the 
dust formation, and (iii) the depletion in the innermost and 
outermost regions. Finally, in Sect. |5] we discuss the results 
and draw some general conclusions. 



2. The chemical model 

The model we are going to develop st ems from the original 
infall model proposed long ago bv (|Chiosil 19801). in the 
lates t multi-ring formulatio n bv iPortinari fc Chiosil (2000) 
and IPortinari et al.l (|2004) for disk galaxies. This latter 
includes the radial flows of matter and the galactic bar 
that were introduced to reproduce the radial gradients 
and the density profile. The adopted stellar yields are the 
original ones by IPortinari et al.l (|l998h in its latest version 
(Portinari 2010 - private communication). The model is 
able to describe the evolution of the abundances of the 
different refractory elements trapped into dust, prop- 
erly simulating th e process of i njectio n of the star-dust 
Q into the ISM (|Piovan et al.l 1201 lal) and dust accre- 
tion/destruction processes whose net balance determines 
the amount of dust in the ISM. The model traces the 
evolution of the abundance of some typical grain families 
supposed to represent the dust in the I SM and single of 
elements trapped into dust ()Piovan et al.ll2011bh . 



The Disk of the MW is subdivided in N concentric circu- 
lar rings where r is the galacto-centric distance. Each ring 

or shell is identified by the mid radius ru with k = 1, , N. 

In most cases, radial flows of interstellar gas and dust are 
neglected, so that each ring / shell evolves independently 
from the others, but in our model the exchange of matter 
between contiguous shell is allowed. The Disk is described 
by the surface mass density as a function of the radial coor- 
dinate r and time t: a (r^, t). Depending on the case, a can 
refer to the ISM (cr M ), in turn split into dust or gas (a D 
or a G respectively), to the stars (<r*) or to the total mass 
(simply a). At every radius r^, the surface mass density is 
supposed to slowly grow by infall of either primordial or 
already enriched gas and to fetch at the present age to the 
mass density profile across the Galactic Disk for which an 
exponential profile is best suited to represent the surface 
mass density distribution: a (r^, to) cc cxp (— r^/Vd), where 
rd is the scale radius of the Galactic Disk, that is typically 
estimated of the order of 4 — 5 Kpc. Since the final density 
profile is a priori known, one may normalize to it the cur- 
rent total surface mass density of the ISM. The evolution of 
the generic elemental species i-th in the dust at the radial 
distance r^ from the centre of the galaxy and at the time 
t, is described as: 



1 Thereinafter, by star-dust we refer to the dust directly pro- 
duced and emitted into the ISM, by AGB stars and supernovae 
(SNae) both of type I and type II. 
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where xf = xY ( r k,t). The first term at the r.h.s. of eqn. 
$5§ is the depletion of dust because of the star formation 
that consumes both gas and dust (assumed uniformly 
mixed in the ISM). The second term is the contribution 
by stellar winds from low mass stars to the enrichment of 
the i-th component of the dust. With respect to standard 
equations for the sole gas component, the condensation 
coefficients 5™ t determine the fraction of material in 
stellar winds that goes into dust with respect to that 
in gas (local condensation of dust in the stellar outflow 
of low-intermediate mass stars). The third term is the 
contribution by stars not belonging to binary systems 
and not going into type II SNae (the same coefficients 
S^ i are used). The fourth term is the contribution by 
stars not belonging to binary systems, but going into 
type II SNae. The condensation efficiency in the ejecta 
of type II SNae are named as S^. The fifth term is the 
contribution of massive stars going into type II SNae. The 
sixth and seventh term represent the contribution by the 
primary star of a binary system, distinguishing between 
those becoming type II SNae from those failing this stage 
and using in each situation the correct coefficients. The 
eighth term is the contribution of type la SNae, where 
the condensation coefficients are named as 8 ci to describe 
the mass fraction of the ejecta going into dust. For all 
th e 5 c .j coefficie n ts we adopt the prescriptions presented 
in iPiovan et al.l (|201Iaf) . The last four terms describe: 
(1) the outflow of dust due to galactic winds (in the 
case of disk galaxies this term can be set to zero); (2) 
the radial flows of matter between contiguous shells; (3) 
the accretion term describing the accretion of grain onto 
bigger particles in cold clouds; (4) the destruction term 
taking into account the effect of the shocks of SNae on 
grains, obviously giving a negative contribution. The infall 
term in the case of dust can be neglected because we can 
assume that the primordial material entering the galaxy is 
made by gas only without a solid dust component mixed 
to it. Equations similar to eqns. ((3J) must be written for 



the ISM and gas. They are identified by the suffix j\4 
(I SM) and G (gas). T hese equations are given in detail 
in lPiovan et al.l (|2011b[ ) and are shortly summarized below. 

Indicating the contribution to the yields by stellar winds 
and type la and II SNae with the symbols Wi,D { r k,t), (and 
similar expressions for gas and ISM as a whole) and ne- 
glecting the outflow term, eqns. (J3J, the equations for the 
the gas and ISM are 



dt 



Mi(r k ,t) = -Xi(rk,t)il>(r kl t) + +Wi,M(r k ,t) 



dt 



Mi(r k ,t) 



rf 
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(4) 



It is soon evident that the dust creation/destruction and 
the radial flows make the system of diff erential equations 
more c omplicated than the original one bv lTalbot fc Arnetn 
(|1975t ) for a one-zone closed-box model. As the ISM is given 
by the sum of gas and dust, only two of these equations are 
required, furthermore Eqn. in the case of radial flows eqn. 
@ can be used only if the gas and dust flow with the same 
velocity. ^^^_^^^^ ^_^^ 

In the model by IPiovan et al.l (|2011b|) , many prescrip- 
tions/values for key quantities are explored. In brief we con- 
sider (i) nine cases for the IMF and five cases for the SFR, 
and for each SFR different specific efficiencies can be used; 
(ii) two possible models of accretion of dust grains in the 
cold regions of the ISM; (ii) different cases of star-dust for- 
mation in AGB stars, type II and type la SNae; (4) the age 
at which the effects induced by the Galactic Bar are taken 
into account. The presence of the Bar is mimicked by a 
change in the pattern of the velocities for the m atter flowing 
from one shell to another (jPortinari fc Chios i 2000). The 
main motivation for the Bar is to reproduce the molecular 
ring at 4-5 Kpc in the gas distribution across the Galactic 
Disk. The velocity pattern is taken from lPortinari fc Chiosil 
(2000) and is suitably chosen for every SFR and age tsar at 
which the e ffect of the Bar is supposed to begin (the reader 
should refer iPortinari fc Chiosil (|2000l ) for all the details). 
In this paper, we keep constant the prescription for t he ve- 
locity patter (Bar) assigned by IPortinari fc Chiosil ((2000) 
to each case of SFR and age tsar- 

Table [T] summarizes the many combinations of the pa- 
rameters and associated models that are possible. For a 
thorough discussion of the effects induced by changing one 
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Table 1. Parameters of the models. Column (1) is the parameter number, column (2) the associated physical quantity, and column 
(3) the sourc e and t he italic symbols are the identification code we have adopted. See the text for some further information and 
I Piovan et al] (|2011bT ) for a detailed description of each choice. 



Parameter Source and identification label 



1 IMF 

2 SFR law 

3 xmc model 

4 Acer, model 

5 SNae la model 

6 SNae II model 

7 AGB model 

8 Galactic Bar 

9 Efficiency SFR 



Salpeter (A), Larson (£>), Kennicutt (C) Kroupa orig. (T>), 

Chabrier (£), Arimoto (J), Kroupa 2007 (g), Scalo {U), Larson SoNe (X) 

Constant SFR (A), Schmidt (B), Talbot & Arnett (C), Dopita & Ryder (£>), Wyse & Silk {£) 

Artificial Neura l Netw orks model (.4), Constan t xmc as in t he Solar Neigh. (B) 

Modified iDwekl ll 19981 ) and iCalura et al] 1120081) (A): adapted IZhukovska et al. I ((20081) model (B) 

Dust injection adapted from: Dwck (1998), Calura et al. (2008]) ( A), IZhukovska et al] (|2008f > (B) 

Dust injection adapted fro m: [Dwck (1991) (A) . IZhukovska et al] (|2008f ) (B), 

iNozawa et aD (|2003l . [2006I . 1200711 (C ) 

Dust injection adapted from: IDwekl l| 19981 ) ( A) . iFerrarotti fc Gaill (|2006h (B) 
No onset (.4), onset at Ig — 4 Gyr (B), onset at to — 1 Gyr (C) 
Low efficiency (A), medium efficiency (B) , high efficiency (C) 



at a time the various parameters on the formation, destruc- 
tion and evolution of (the Sola r Neighborho o d, in p articu- 
lar) the reader should refer to iPiovan et al] (|2011bl) . High 
number of parameters implies a high number of possible 
models. So a code must be found to establish a easy to 
read correspondence between the particular combination of 
parameters and the associated model. Each model is iden- 
tified by a string of nine letters (the number of param- 
eters) in italic face whose position in the string and al- 
phabetic order corresponds to a particular parameter and 
choice for it. The code is presented in Table [TJ The se- 
quence should be read from top to bottom: for instance, the 
string VBAABABAB corresponds to Kroupa 1998 IMF, 
Schmidt SFR, ANN mod el for the fra ction of molecular 
clouds in the ISM xmc IDwekl ([1998) accretion model, 



Zhuk ovska et al. (120081) r ecipe for the yields of dust from 
type la SNae, Dwek (1998) condensation efficiencies for type 
II SNas, IFerrarotti fc Gail (2006) condensation efficiencies 
for AGB stars, no effect of the Bar on the inner regions, 
and high efficiency v of the SFR. Radial flows are always 
be included by default. We do not explore here the effects of 
diffe rent IMFs a nd/SFR that have been already examined 
bv IPiovan et all (|2011bD . 

We focus here on the radial behaviour of dust formation 
and evolution and seek to answer the following questions. 
What are the dominant types of dust at different ages of 
the Galaxy and radial distances? Are the predicted radial 
gradient of element abundances in agreement with the ob- 
servations in such a way to allow a correct calculation of 
dust formation on consistent initial conditions? What are 
the differences between the inner and outer regions of the 
Galaxy? To highlight these and other related issues we set 
up the so-called "reference model" in which we vary the key 
parameters. This reference model is identified by the string 
QVABBCBBB and it includes the most refined recipes for 
the production of dust by stars and the accretion of dust 
in the ISM. It is worth noticing here the combination of 
the parameters we have chosen does not necessarily yield 
the model best fitting the observations across the Galactic 
Disk, but it only meant to provide the gross features on 
the physical processes for dust formation across the Disk 
and clues on the general behaviour of the system. Clearly, 
the model must be in satisfactory agreement with the ob- 



servational data, even though it does not perfectly match 

it. 



3. Modelling radial flows 

A key term in the description of the properties of gas and 
dust along the Disk is the one of radial flow in eqn. driving 
the exchange of gas/dust between contiguous shells. The 
ra dial flows are describ e d acc ording to the model developed 
by iPortinari fc Chios! (|2000t ) however extended to a two 
components fluid made of gas and dust each of which in 
principle moving with its own velocity. Recalling that each 
circular ring /shell is identified by the mid-point radius r k 
and the inner and outer borders by the radii r k _i and r k , i, 
the material can flow across the inner and outer borders 



with velocity v k _i and 



'k+i 



respectively. Flow velocities 



are taken positive outward and negative inward. Let denote 
with Fd (r) the radial flow of dust and Fq (r) that of gas, so 
that the total flux of ISM is F M (r) = F D (r) + F G (r). The 
motion of this two-components fluid alter the gas surface 
density in the fc-th shell. The most convenient description 
of it is as follows: If there is no velocity drift between gas 
and dust, as we will later assume, it is more convenient to 
follow the ISM and the dust, whereas for different gas and 
dust velocities it is more convenient to follow dust and gas 
as two independent fluids. Therefore we have: 



da c r (k) 



dt 



1 



rf 



7T I 7 



k+4 



AF C (r fc ) (5) 



where C = A4,G,D (C is for component) means that we 
have the same formal expression for ISM (A4), gas (G) and 
dust (£>). The term A is F c (r k ) = F c (r k+ i) - F c (r k _i). 

The flux across the external border of the ring/shell r k+ i 
is: 



F c (r k+ i) =2irr 



k+ 



h V Uh 



x{v C k+h )<? C K + 



+ x(-v e k+h )4+i 



(6) 



where x( x ) ^ s the step function: x( x ) = 1 or for 
x > or < 0, respectively. Eqns. (JU are a sort of "upwind 
approximation" for the advection term to be included in 
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the model equations, describing either inflow or outflow de- 
pending on the sign of the velocity v9, 1 . Clearly, since 

there are two components, gas and dust, it may happen 
that v? i 7^ v? i . The simplest case is for v G — °' D — 



k+4 



k+ 



k+4 



k+4 



c 



i , i.e. when gas and dust are coupled and move along 

2 

the same direction. Analogous expression holds for the flux 
Fc(r k _i) across the inner border of the ring/shell. 

Let's take the inner edge r k _ i at the midpoint between 
Tk—i and r k , and similarly for the outer edge r k+ i. Writing 
Eqn. ([5]) separately for each chemical species i and as func- 
tion of Gj's, we may write the radial flow term in eqn. [3] 
and similar ones for other components as follows: 



be kno wn in advance. It is know n indeed only a posteri- 
ori (see iPortinari fc Ch iosi 2000, for more details). At the 
end of each simulation we need to check how much radial 
flows have altered the actual density profile a{r k ,tc) with 
respect to the pure accretion profile ovi(rfc). In any case, if 
the flow speeds are small (of the order of v ~ 1 km sec -1 ), 
the two profiles are expected to be close each other. 

3.1. Boundary conditions 

Eqns. need to be modified for the innermost and the 
outermost shells, since for these regions at shells fc— 1 or 
fc+1 are not defined. 



cit 



Ci(r k ,t) 



= a£e i (r fe _i,t) - /3%Ci(r k ,t) + 

rf (7) 

+ lkC t (r k+ i,t) 
where again C = Ai/G/D and: 



Ofc = 



rt + 



r k-i+ r k+i 



, c v c r k -! +r k 
2 2 ^fc+i — r k _i 



0A(k-l) 



OAk 



0k = - ■ 



rk + 



r k-i+ r k+i 



X{-v k _i)v k _ 



r k -i +r k 
' r k +i - r k -i 



-X(v k+ i)v^, 



rk +r k +i 

x 

2 r k +i — rfe_i 
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PAk 
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(8) 



The terms at the r.h.s. of the eqn. ([7]) describe the con- 
tribution of the 3 contiguous shells to driving the mat- 
ter flows: the first term represents the gas being gained 
by shell fc from the shell fc— 1, the second term is the gas 
lost by the shell fc towards the shells fc-1 and fc+1, and 
the third term is the gas acquired by the shell k from 
the shell fc+1 (jPortinari fc Chiosil [2000h . The quantities 
°\4(fc-i)i o~Ak, and o~A(k+V) are compact notations for the 
present-day surface mass density profile a(rj,to) = <JA(rj) 
where j = fc — 1, fc, fc + 1, and A stands for accreted matter. 
The coefficients given by Eqns. © are all > and depend 
only on the shell fc, not on the chemical species i considered. 
They will be different for gas and dust if a drift between gas 
and dust is considered. If there is no drift, then the M. equa- 
tion can be used together with one between dust and gas 
equations. If the velocity pattern is constant in time, a k , 
f3 k and 7^ are also constant in time. There is an important 
point to note concerning the present-day surface mass den- 
sity profile. In the case of static models (no radial flows) the 
present-day surface mass density is known a a priori and is 
determined by the mass accretion law one has assumed, i.e. 
a{r k ,tc) = o~A(rk) = OAk- In other words, in static mod- 
els the radial profile for accretion can be directly chosen 
so as to m atch the observed pre s ent-day surface density i n 
the Disk (|Portinari et alJll998l IPortinari fc Chiosil fl999h . 
The inclusion of radial gas flows changes the expected fi- 
nal density profile (as it is now set up by infall and radial 
flows). Therefore, a(r k ,tc) ^ CA^k) and a(r,tc) cannot 



3.1.1. The innermost shell 

Our model is limited to the Disk of the MW and can- 
not follow the innermost regions of the Galaxy where the 
Bulge dominates. A correct description of this region would 
require a more complicated model: a spherical compo- 
nent simulating the Bulge combined with the Disk each 
of which with its own mass distribution and formation his- 
tory. Therefore we truncate the Disk to a innermost bound- 
ary (ring) located where the Bulge stars dominating, i.e. at 
about ri = 2 ~ 2.5 kpc. At this innermost boundary, the 
first shell is taken symmetric with respect to r% and its 
innermost radius is 



•i=(3n-ra)/2 



(9) 



At this layer we impose that always v^< 2 < 0, therefore 
fc = 1, and Eqn. ([7]) become 



dt 



Ci(n,t) 



-tfCi(n,t) + ^Ci(r 2 ,t) 



(10) 



rf 
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3.1.2. Boundary conditions at the Disk edge 

The Galactic Disk is conceived as made of a inner region, 
where star formation and chemical enrichment occurs and 
stars gas and dust exist, and a external region in which 
only gas is present, no stars can be formed and no chem- 
ical enrichment can take place. This picture is supported 
both by observational data showing that in external spirals 
HI disks are observed to extend much beyond the optical 
disk and theoretical considerati ons about stability prevent - 
ing SF beyond a certain radius (|Toomre| |l964: Qu irlall972l) . 
The part of the Disk we are interested in is the first one, 
which may extend up to about 20-22 kpc. Therefore the 
mid radius of the most external shell (fc = N) can be lo- 
cated at rw = 20 — 22 kpc. all the regions belonging to 
the Disk external to this limit can be considered as a large 
reservoir of gas with primordial composition from which gas 
can flow inside. At the present time when the gravitational 
settling of the proto-galactic cloud is over, the radial flow of 
gas from this external region of Disk can be even more effi- 
cient than the classical infall mechanism. If no star forma- 
tion occurs, the evolution of the gas surface density in the 
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outer disk Vr > rjy + i/2 can be described by (jLacev fe Falll 
ll985HPortinari fc Chiosill200fti : 



d[G(r,t)} 



1 d 



= A(r)e^ -LJL[rv.G(r,t)] (12) 

at r or 

Without star formation, there is no chemical evolution 
and the pattern of abundances always remain the primor- 
dial one (Xi^nf). We assume that this primordial material 
is also dust-free, that is D(r,t) = for r > r 



N+4 



As 

a consequence of it, the normalized surface density of the 
outer Disk is given only by M(r 7 t) — G{r, t) for r > r N+ 1 

and this component is the only one to be considered in Eqn. 
(fT2"|) . Other simplifying conditions for th e the most external 
part o f the Disk have been discussed by iPortinari fc Chiosl 
J|2000D . They all apply for r > r N+ i. In brief, the infall 
time-scale is uniform r(r) = r(r]\[); the inflow velocity is 



uniform and constant v(r, t) 



J N+i 



and the infall profile 



is flat A{r) = A ext in agreement with the observed gas discs 
in spirals. With these assumptions, Eqn. (fT2"j) becomes: 



d[G(r,t)} d[G(r,t)] 



Ae 



v 

— a 
r 



(13) 



dt dr 

where r = t(tn), v = v N+ i, and A = A ext to simplify the 
notat ion. Eqn. (|13p has a stra ightforward analytical solu- 
tion (|Portinari fc Chiosil[2000l) : 



a(r,t) = At x 
l-e~ 



- (t - Trf, 



(14) 



where T r f > is the time when radial inflows are assumed 
to start. Eqn. (fT4"|) is our boundary condition at the outer- 
most edge. 

Notice that relation (|14j) is the solution of Eqn. (fl~3|) 
in the idealized case of an infinite, flat gas layer extend- 
ing boundless to any r > r^ (see also Appendix B). Of 
course, this does not correspond to gaseous disks surround- 
ing real spirals. However, as we are considering only slow 
inflow velocities (v ~ 1 km sec -1 ), with typical values of 
tn = 20 — 22 kpc and to = 12.8 Gyr, the gas actually drift- 
ing into the inner disk shells will be just the gas originally 
accreted from regions within r ~ 35 — 40 kpc. Therefore, 
the boundary condition (|14p remains valid as long as the 
gas layer extends out to ~ 35 — 40 kpc, a very plausible 
assumption since observed gaseous disks extend over a few 
tens or even ~ 100 kpc. 

3.1.3. The outermost shell 

We take a reference external radius r ext > r^ in the outer 
disc where the (total and gas) surface density u{r ext ,t) = 
Cext(i) is given by the boundary condition (TT4")) : typically, 
rpj = 22 kpc and r ext ~ 24 kpc. We take the outer edge of 
the shell at the midpoint: 

?> + i = (r N + r ext )/2 (15) 

The radial flow term for the N-th shell is: 



—C^rN^t) 
dt 



rf 



=a N Ci(r N -i,t) - P c N Ci{r N ,t) + 

U f(t) (16) 



where: 
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(17) 



rjv-i+2rjv+r ext r ext — rjv-i 



for M and G, while for dust wf (t) = 0. 



3.2. The numerical solution 

Using Eqns. 0, PH|) and ([T^]). neglecting for the moment 
the terms for dust accretion and destruction, the basic set 
of Eqns. ([2]), ([3]) and (HJ can be written as: 

( jMrut) = tfi (t)Ci(n,t) +7i Ci(ra,t) + Wi,c(n,t) 

d 



irCi(r k ,t) =aiC i (r k -i,t)+'9i{t)Ci{r kl t) + 

+ liC t {r k ,t) + W l ,c{r k ,t) 
k = 2,...,N- 1 

-^Ci(r N , t) = a%Ci(r N -i,t) + $ c N (t) Ci(r N ,t) + 

+ W it e{rN,t) + cjf(t) 

where we have introduced: 



(18) 



m) 



(v C (r k ,t) 



* 



Pi) 



<0 



V ( r k,t) = T {r k ,t) 



(19) 



(20) 



The terms Wi c(rk,t) are defined as in Eqns. (O, ([3]) and 
dU). We refer to [Portinari et all ()1998t ) for further details 
on the quantity r\. Neglecting, for the time being, that ry's 
and Wi's contain the Cj's themselves, we are dealing with a 
linear, first order, non homogeneous system of differential 
equations with non constant coefficients, of the kind: 



dt 



= A(t)d(t) + Wi(t) 



(21) 



There is a system like (f2i"j) for each chemical species i, but 
the matrix of the coefficients A(t) is independent of i. The 
way by which this system is solved and then dust evolution 
is calculated depends on the specific hypotheses made for 
gas and dust. 

(1) Radial flows are neglected but dust is included. Since 
the radial flows are not started from the beginning of the 
evolution of the galaxy, this case usually happens in the 
early stages of the evolution of the MW and it lasts some 
Gyrs (jPortinari fc Chiosl 12000') . In this case the evolution 
of the M. and D components should be followed. The 
classical static system for the ISM in which the different 
rings evolve independently ([2]) (without flowing term) is 
solved according to th e classical Talbot & Arnett technique 
(jPortinari et al.lll998l) . Once this system is solved the accre- 
tion/destruction terms for dust can be taken into account 
(see Piovan et al. (2011b)) is solved to find the evolution of 
the dust. The presence of the accretion/destruction terms 
does no longer allow to use the standard resolving method 
and the new system of ordinary differential equations must 
be integrated with a different technique. We adopt a com- 
bination of the ODEINT r outine (4-th order Runge-Kutta) 
with controlled step-size (|Press et al.l 119921 ) and the pow- 
erful DOP853 Runge-Kutta routine (method of Dormand 
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& Prince) of ord er 8 with a 5th orde r estimator and a 3rd 
order correction (jHairer et al.l 1201 Pa lo). 

(2) Radial flows and dust are both included. The equations 
to be solved depend on whether or not some velocity drift 
exists between dust and gas. The simplest and most widely 
adopted hypothesis is that not only dust is homogeneously 
mixed to the gas, but both radially flow at the same veloc- 
ity. In this case, the two systems of equations to be solved 
are © and ©. The system © is reduced to the form (TT8"|) 
by including the radial flows and solved accordi ng to the 
technique developed bv lPortinari fc Chiosil ((2000) to whom 
the reader should refer for all details. Once the evolution of 
the ISM is computed, Eqn. ([2]) for dust is setup by means 
of accretion/destruction terms (|Piovan et al.l 1201 lbft . The 
system of ordinary differential equations is now solved with 
the ODEINT and DOP853 integrators. The reason to cal- 
culate first the evolution for the ISM is that the equations 
for the accretion/destruction processes regulating the dust 
evolution include the total amount of interstellar matter as 
an ingredient. 

(3) Velocity drift between the ISM components. In such 
a case, the global system of differential equations, now in- 
cluding gas and dust together, made by Eqns. dH) and ([3]), 
with the radial flow term as in Eqn. ([7]) and the accre- 
tion/destruction terms for dust must be solved. In such a 
case it is not possible to obtain a priori the evolution of 
the ISM. Again the ODEINT and DOP853 integrators are 
used to evolve the system as a function of time. 

Some general considerations can be made: rather small 
time-steps are needed to keep stable the numerical model; 
the required time-steps get smaller and smaller at increas- 
ing the flow velocities and making thinner the shells. The 
integrators of ordinary differential equations self-regulate 
the internal time-step according to the required precision 
once the external time-step is chosen. To describe the gas 
flows in a disk with an exponential density profile, the shells 
are equi-spaced in the logarithmic, rather than linear, scale 
(so that they roughly have the same mass, rather than the 
same width). We model the Galactic Disk using 33 shells 
from 2 to 21 kpc, equally spaced in the logarithmic scale; 
their width ranges from ~ 0.15 kpc for the inner shells to 
~ I — 1.5 kpc for the outermost ones. With such a grid spac- 
ing, and velocities up to ~ 1 km sec , suitable time-steps 
are of 10 -4 Gyr. This means that roughly 1.5 x 10 5 time- 
steps, times 35 shells, are needed to complete each model, 
which would translate in excessive computational times. 
This drawback is avoided by separatin g the time-scales in 
the c ode following the suggestion by (jPortinari fc Chiosil 
l2000ft . 

4. The Radial variation of the fraction of MCs 

The accretion of dust in the MW Disk can be simulated 
in our model by means of two possible choices: (1) a sim- 
ple model that adopts typical accretion timescales to sim- 
ply determine the amount of each element embedded into 
dust; (2) a more refined model with variable timescales that 
follows the accretion of dust grains in some typical grain 
families inside the cold regions of the Disk. This latter de- 
scription is the one included in our reference model and 
is adopted here to examine the radial behaviour of dust 
formation and evolution. Since in our one-phase chemical 
model no description is available for a multi-phase ISM, it 
is unavoidable to estimate from time to time (possibly re- 
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Fig. 1. Upper panel: present-day radial profiles across 
the MW Disk for molecular hydrogen an 2 (r,tc), neu- 
tral hydrogen <jhi (r, to), and total amount of gas 
a (r,ta) obtained taking into account h elium correction 
(jBoissier fc Prantzosl[l99a I Yin et alJl2009h . Lower panel: 
present-day radial profile of star formation across the 
MW Disk according to the data by iGuibert etall d 1978ft 
squar es), iGuesten fc Mezgen ( 11982') (circles), iLvne et all 
1985) (hexagons). iRanal (|l99ll ) (stars). Superposed to the 
data are three theoretical SF laws (Schmidt - continuous 
line, Wyse & Silk - dotted line, Dopita & Ryder - dashed 
line) and an exponential best fit, all referring to the present 
age t G of the MW. 



latcd to the physical properties of the system) the amount 
mass in cold molecular clouds (MCs) in each ring. The par- 
tition between MCs and diffuse ISM (the remaining part of 
the ISM) crucially enters the accretion term in eqns. [21 El 
and E] and drives the evolution of the dust budget. In this 
formulation the fraction, xmc, of ISM locked up in MCs is 
a parameter. For this reason we have included two possibil- 
ities for \mc- In the first one, f r om th e present data for the 
Solar Vicinity [Zhukovska et al.l ([2008D estimate xmc = 0.2, 



and assume that it has remained constant from the for- 
mation of the MW to the present. Therefore, a constant 
value is adopted during the entire evolution of the Solar 
Neighbourhood. The second description for xmc is pre- 
sented below and it stands on the data for the MW Disk. 
In practice, it correlates the fraction of MCs to t he SFR 
and the local gas density. In iPiovan et al.l (|2011bft we ob- 
served as both recipes produced similar results for the Solar 
Vicinity, but in the case of a constant fraction xmc the 
amount of dust could be underestimated and overestimated 
in the innermost and outermost regions of the MW disk, 
respectively. If compared to other more physically sounded 
prescriptions, a constant xmc would give rise to a flatter 
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"gradient in dust". We expect however that the local phys- 
ical conditions affect the accretion process and the fraction 
of MCs varies with time and space. This strongly favors the 
second prescription for a varying fraction of MCs. 

Fig. Q] shows the radial dependence across the MW 
Disk of the current surface density profiles of H2, HI 
and total gas (upper panel) of the MW and the current 
SFR normalized to the solar value SFR/SFRq (lower 
panel) indicated by the observational da ta. Superpose d to 
the da ta we show some typical SF laws iPiovan et al.1 (See 
l2011bL for more details about these analytical laws). It is 
worth recalling that the total gas sur face density includes 
also the contr i butio n by helium (|Boissier fc Prantzosl 
119991: lYinet al.l [20091 ). With the aid of this, we look at 
the correlations between SFR/SFR Q , <j G (?*&,£) (one of 
quantities followed by the chemical model), and an 2 i r k,t) 
(related to the parameter xmc)- 

To this aim, first we derive the best fit for the three 
quantities: the results (obtained using the MATLAB Curve 
Fitting Toolbox) are indicated by the solid lines in the top 
panel and the thick line in the bottom panel of Fig. [1] 
Second, using the best fits we derive the 3D correlation be- 
tween the H2 mass, the star formation rate, and total gas 
mass. They are shown in the three panels of Fig. [2] To repro- 
duce the 3D data and extrapolate the results to unknown 
regions of the parameter space, we use the Artificial Neural 
Networks technique (hereafter ANNs) modelled from bio- 
logical nervous systems: ANNs are like a network of neu- 
rons linked together by synapses (numerical weights). The 
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(2008) (Se08). iMagrini et al.l l|2009l ) (Ma09); Cepheid stars: 
An0 2a, An02b, An 2c, An 02d, Lu03, Ko05, Lu06, Y0O6, Le07 
and ILemasle et al.l l|2008t ) (Le08 ). Carbon. O and B stars 
(Gu98, SmOl, Da04); HII regions: lEsteban et al.l ((20051 ) (Es05); 
Cepheid stars: An02a, A n02b, An02c, An02d, Lu03, Ko05 and 
lAndrievskv et al.1 (2004) (An04). Silicon. O and B stars Gu98, 
SmOl, Da04; Field Red Giants: Ca05; Open Clusters: RoOO, 
Ca04, Ca07, Yo05, Se08; Cepheid stars: An02a, An02b, An02c, 
An02d, Lu03, An04, Ko05, Lu06, Y0O6, Le07. Sulfur. Cepheid 
stars: An02a, An02b, An02c, An02d, L u03, An04, Ko05, Lu06; O 
and B stars: S mOl, Da04; HII r egions Martin-H ernandez et al.l 
(p002T ) (Ma02), iRudolph et all ((20071) (Ru07); Open Clusters: 
Ca04, Ca07, Yo05, Se08; Field Red Giants: Ca05; Cepheid 
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Y0O6, Le07. Calcium. Open Clusters: Ca04, Ca07, Yo05, Se08; 
Field Red Giants: Ca05; Cepheid stars: An02a, An02b, An02c, 
An02 d, Lu03, An04, Ko05, Lu 06, Y0O6, Le07. Oxygen: HII re- 
gions (Dehan^engltal. (2000) (DeOO), Es05, Ru07; O and B 
stars: ISmartt fc RoUestonl IJ1997D (Sm97), Da04, SmOl; Field 
Red Giants: Ca05; Open Clusters: RoOO, Ca04, Yo05; Cepheid 
stars: An02a, An02b, An02c, An02d, Lu03, An04, Ko05, Le07. 
Nitrogen: HII regions: Ru07; O and B stars: Gu98, Da04, SmOl; 
Open Clusters: RoOO; Cepheid stars: An02b, An02c, Lu03, Ko05. 



network in usage here is taken from iGrassi et al.1 ([2011) 
to whom the reader should refer for all details. It con- 
tains two layers of hidden neurons whose number is rih — 5 
and n g = 25, respectively and has two input neurons, i.e. 
(SFR/SFR Q and and a G (rk,ta)) (plus the bias input that 
is always set to a constant value) and one output neuron, 
i.e. cth 2 (^fej^G)- The evolution of the RMS is monotonic 
and decreases to approximately RMS = 10~ 3 after 10 6 it- 
erations. In the left panel of Fig. [5] we show the fit found 
by the ANNs (solid line) compared with the data (crosses). 
We suppose now that the correlation between the three pa- 
rameters established from the present-day data remains the 
same as a function of the Galaxy age and therefore we drop 
the dependence on time. Furthermore, adopting the same 
ANN weights, we extrapolate the ANN predictions to the 
whole space of the parameters SFR/SFRq, a G (r,t) and 
(T{j 2 (r,(). The contours and isolines describing the extrap- 
olated regions are shown in left and middle panels of Fig. 

El 

Some comments are worth being made here. First, we 
have no warranty that the behaviour in the extrapolated 
region is so regular as predicted by the ANNs. In the huge 
gaps uncovered by the experimental data for the MW, we 
could for example have valleys and peaks so that to cast 
light on this problem a complete coverage of the parameter 
space with observational data for spiral galaxies in general 
would be required. However, this goes beyond the aims of 
this study. It suffices here to check that extrapolating the 
MW data does not lead to un-physical results. The predic- 
tions seem to be reasonable: for low SFR or low gas amount, 
the H2 mass is low. For high SFR and gas mass we have a 
reasonable amount of H2 . Even if we look only at the region 
where the SFR is high but the total amount of gas is low, 
the prediction for H2 is also reasonable. Even if the corre- 
lation between H2 and SFR is not clear and somewhat in 
contradictions with the expe ctations, stars ar e seen to form 
from the cold molecular gas (jKennicuttJ 1 1 9981 ) and a signif- 
icant amount of H2 is physically acceptable when SFR is 
significant. 

Second, the region covered by the ANN is bounded 
by the limits [(SFR/SFR ) mm , (SFR/SFRq ) ma ,] and 
[a G (r, t) min , a G (r, t) max \. This may be a problem in par- 
ticular with high SFRs, because the limits for the ANN 
have been set basing upon data for the SFR in the inner- 
most regions of the Disk where the uncertainty is high. To 
cope with this point of embarrassment, when the SFR is 
higher than about 3 x SFRq, we replace the predictions 
based on the ANNs with a simpler fit based on ajj 2 ( r , to) 
versus SFR data and we keep the fraction xmc in the in- 
terval 0.4 — 0.5. All of this is shown in the right panel of 
Fig. El 



5. Radial gradients in element abundances 

Studying the properties of dust in the SoNe is by far easier 
than across the Disk of the MW: indeed the data to our 
disposal for the Disk are much fewer and of poorer quality 
than those for the SoNe. For instance, we cannot estimate 
the depletion of elements at different distances from the 
SoNe toward the Galactic Center and/or the outermost 
regions of the Disk. However, we can compare the radial 
gradients in element abundances, as indicated by Cepheid 
stars, Open Clusters, HII regions and 0,B stars. In the 
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Fig. 2. Left Panel: 3D plot of <JH 2 ( r >tG) a s a function of SFR/SFR and cr G (r, £g) for the present-day radial profiles 
across the MW Disk (crosses). The solid line superposed to the data is the fit from the ANNs. Contour isolines obtained 
by extrapolating the ANNs weights to the regions uncovered by the data are also shown. Middle Panel: The same as 
in the Left Panel but showing only the 2D contour levels of <7H 2 { r , to) hi the SFR/SFR© vs. a G (r,tc) plane. Right 
Panel: 2D plot with the radial SFR/SFR© vs cr// 2 (r, tc)/<y G {r 1 tQ) relationship for the MW (triangles) and a two-piece 
fit of the data. The region of interest is the one for SFR higher than solar. 
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Fig. 3. Radial surface density profile of the gas in the MW for twelve models taken from iPiovan et all (|2011bh . The 
models differ for the parameters: IMF (Larson SoNe, Salpeter, Kroupa multi power law), infall timescale r (3, 6 and 9 
Gyr), the efficiency v (0.3 and 0. 7). Not all th e combinations have been represented, but only the most significant ones. 
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the errors in the ob s ervati onal estimates. Data for the SoNe gas mass density are from iDickevI (|1993| ): lDamel (|1993l ) and 
lOlling k Merrifieldl (f200lh 



case of a satisfactory agreement between the observations 
and theoretical predictions for the element abundances, we 
may argue with some confidence that the predictions for 
radial variation of the dust properties are reasonably good. 
In relation to this see the case of the relative abundances 
of Mg and Si and how they would affect the formation 
of silicates in the dust and depletion of thes e elements as 
thoroughly discussed by IPiovan et al.l (|2011bh . 

In the following, we examine the radial variations of the 
elemental abundances to check if our multi-shells model, 
in presence of radial flows of gas and inner bar, is able to 
reproduce in a satisfactory way at least the main radial 
properties of the MW. As already said we are not looking 
for the best model in absolute sense but only the a general 
agreement with observational data so that our predictions 



for dust formation and evolution stand on solid ground. 

Two important observational constraints must be 
satisfied and reproduced by the models, namely: (i) the 
radial distribution of the gas and (ii) the radial gradients 
in those elements that are involved in the dust formation. 
In Fig. [31 starting from the reference model introduced 
in Sect. [2 we compare the density profiles for models 
calculated with different IMF, efficiency of the star 
formation v and the infall timescale r. Twelve models 
are shown. This sam ple is taken f r om the wider grid of 
models calculated by IPiovan et al.l (|2011b| ) to reproduce 
the depletion of the elements observed in the SoNe q The 



2 The whole sample contains model for four vt (t = 1, t = 3, 
t — 6 and r = 9 Gyr), three IMFs (Salpeter, Kroupa multi- 
power law and Larson adapted to the SN), four values of the 
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twelve models ar e shown in Fig. El The y are with original 
stellar yields by iPortinari et al.l (|1998) (no correction to 
the Mg abundance is applied). Since we have no data for 
the element depletion as a function of the radial distance, 
the slight under-abundance of Mg is less of a problem here. 
Furthermore, selecting the twelve cases, very short infall 
timescales have been excluded because the models would 
not reproduce the element depletions observed in the 
SoNe. Finally, since we have no hints on the observational 
error associated to the mass density profile cr g (r) as a 
function of the galacto-centric distance, at each distance 
we apply an error derived from the ratio Aa g /a g for the 
SoNe, and assumed to remain constant. Indeed, already 
in the SoNe, there is a significant uncertainty in the 
local gas density, the estimates of which g o from 7-13 
M^pc" 2 (lDickevlll993D. t o ~ 8 (|Damd 119931) up to 13-14 
(|01hng fc Merrifieldll200ih . 

It is soon evident that while some models fairly agree 
with the observational data for a g , others significantly 
deviate from it in particular for the innermost regions: 
models with long infall timescale and low star formation 
efficiency (such as the model with r = 9 Gyr and v = 0.3) 
do not consume enough gas. In literature the problem is 
often solved either introducing a suitable radial dependence 
for t or changing the IMF in such a way that the low-mass 
long-lived stars are favored. However, for the purposes 
of the present study, we prefer to keep both r and IMF 
constant with the galacto-centric distance. Finally, it is 
worth noting the crucial contribution of the bar whose 
presence nicely reproduces the gas density distribution in 
the innermost regions, simulating the bump at about 5 kpc. 

In Fig. 2] we show the gradients across the Galactic Disk 
in the abundance ratios [Mg/H], [Fe/H], [C/H], [Si/H], 
[O/H], [S/H], [Ca/H] and [N/H]. All these elements are also 
involved in dust formation and evolution. For each element 
we group the observational data according to the source of 
the abundances, namely O and B stars, field Red Giants 
(RGs), HII regions, Cepheid variables and open clusters. 
The references for the observational data are listed in 
the footnote 7 . Only models with intermediate values of r 
(1 < t < 9) are compared with the data. Since most likely 
different sources trace the gradients at different galactic 
ages, we plot the gradients as they were at the time the 
Sun was born (about 4.5 Gyr ago, the thin blue and black 
lines) and at the present age (thick red and magenta lines). 
Most objects under examination should fall into the region 
spanned by the gradients during this time interval. In 
general, theoretical predictions and data fairly agree with 
minor discrepancies that can be, in some cases, attributed 
to the yields, such as the slight Mg under-abundanc e 
and/or C and S over-abundance (jPortinari et al.1 [T_998) . 
Other major discrepancies between the theoretical and 
observational abundance gradients that depend on the 
type of source under consideration are expected and can 
be taken to indic ate the pe rforma nce of the models. For 
instance, following iKovtvukh et al.l (|2005|) . in intermediate 
mass stars the C under-abundance (carbon is deficient 



efficiency of the star formation u (y — 0.3, v — 0.7, v = 1.1 
and v — 1.5), and two cases for the stel lar yields of elements, 
the original ones bv lPortinari et all l| 19981 ) and a revised version 
taking into account a correction for a slight under-abundance in 

Mg. 



at about [C/H] = —0.4) and the N over-abundance 
(nitrogen is e nhanc e d at a bout [N/H] = +0.4) h a s been 
predicted by iLuckl (|1978[ ) and iLuck fc Lambert! (|1981l) 
and confirmed in studies of Ce p heids a nd non-variable 
super -giants (|Andrievskv et al.1 Il996t IKovtvukh et al.l 
1996). This is more or less what we observe when the 
evolution of theoretical element abundances C and N along 
the Disk are superposed to the measurements of C and N 
abundances in Cepheids (intermediate mass stars of 3 to 
9 M©). Finally, comparing theory and data, one should 
keep in mind that models do not take star dynamics and 
kinematics into account and the stars always remain in the 
same region where they are born: the radial flows act only 
on gas and dust. 

Radial gradients in dust. Assessed that the radial gradi- 
ents, in the gaseous element abundances and total gaseous 
density, reasonably agree with the observational data, 
at least for some combinations of the input parameters, 
we examine the radial formation and evolution of dust 
in the MW Disk being confident that it is grounded 
on a realistic evolution of the global abundances in the 
ISM. In Fig. \E\ we show for five selected ages, from the 
early stages to the current time, the time evolution of 
the radial abundance of four typical grain families, that 
is silicates, carbonaceous grains, iron grains and generic 
grains formed with Ca, N and S. We grouped these last 
three elements into a single group for the sake of simplicity 
because of their low abundance. The contributions to 
each one of the four groups have been split among the 
three sourc es of dust, that is star-dust from SNae and 
AGB stars (|Piovan et al.l 201 laf) and acc retion in the cold 
regions of the ISM ( Piovan et al.l 1201 lbl ). The amount of 
dust in all the regions continuously grows following the 
enrichment in metals produced by the many generations 
of stars. In the early stages of the evolution, at 0.1 Gyr, 
star-dust from SNae is the main source of dust across the 
whole disk. At 0.5 Gyr the accretion of dust in the ISM 
starts to be significant in the innermost regions and slowly 
becomes the dominating source in the whole Disk, earlier 
for the innermost regions and later for the outermost ones. 
AGB stars contribute mostly with carbon-rich grains in 
the early evolution at low metallicity Z and only when 
a significant metallicity is reached they start enriching 
the ISM in silicates. It is worth noting here that the 
evolutionary scheme across the Disk we have just outlined 
is a consequence of the role played by the star-dust from 
SNae in the early stages. Reducing the yields of dust by 
SNae would obviously change the the above results for the 
early stages. 

In Fig. [S] we present the evolution with time of the 
radial dependence of the logarithm of the dust abundance 
af (rk,t) for all the elements in our list taking part to 
dust formation, namely C, N, O, Mg, Si, S, Ca and Fe. 
Nine ages are represented from the early stages to the 
current time. All the elements more or less follow the same 
evolutionary trend, and eventually reach a typical profile 
with a bump at 4-5 kpc in agreement with the gas density 
profile (as expected). It is interesting to note that while 
for the innermost-central rings of the disk the enrichment 
in dust already reaches a maximum at about 5 Gyr, for 
the outermost regions the dust budget keeps growing 
until the present age. The reason for that can be easily 
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attributed to the delayed contribution by accretion in the 
ISM occurri ng in the oute r regions with low star formation 
(see iPiovan et al.l l2011bl for more details). In brief, as 
accretion of dust grains in the ISM gets important only 
when and if some enrichment in metals has taken place, in 
the outer regions of low star formation this is inhibited by 
the low metal content that can be reached over there. 

It is interesting to examine not only the abundance 
evolution of each element in the dust, but also the relative 
abundance of each element embedded into dust (or the 
relative fraction of each dust grain family) with respect to 
the total dust budget. For this reason in Fig. [7] we show 
the evolution of the radial mass fraction of each element 
embedded into dust o~P (rfe , t) j a D (r^ , t) , normalized to 
the total dust budget, for some of the elements in our list 
and/or some families of grains. In particular, we show 
the fraction of silicates (which means Mg, Si, Fe and O 
embedded into quartz/pyroxenes/olivines by accretion in 
the ISM or star-dust) and carbonaceous grains (this latter 
is nearly coincides with the amount of C, once subtracted 
the amount of C belonging to the silicon carbide). Let us 
examine the behaviour of silicates and carbonaceous grains 
in some detail. At the beginning of the evolution, where 
only the most massive SNae contribute, the mass of dust in 
silicates is a bit larger than in carbon based grains. This 
is in agreement with the recent observatio ns of SN 1987A 
with the Herschel Space Observatory (jMatsuura et al.l 
1201 ll) : the significant emission by a population of cold dust 
grains observed by PACS and SPIRE can be reproduced 
only with a mixture of both carbonaceous grains and 
silicates with these latter more abundant than the C-based 
dust grains. This is indeed what is shown in Fig. [S] at 
the age of 0.1 Gyr, when SNae play the major role. As 
soon as AGB stars start to contribute (at low/sub-solar 
metallicities a 5-6 solar masses has a lifetime of the order of 
0.1 Gyr) the mixture of SNae and low-Z AGB stars slowly 
increases the fraction of carbonaceous grains and reduces 
that of silicates. Indeed, low-Z AGB stars are easily C-star 
and mainly inject carbon in the ISM. At ages in the range 
0.5-1 Gyr, at least in the innermost regions, the formation 
of dust grains in the cold ISM starts to be a very efficient 
process and strongly influences the total budget. Our 
silicates are made by Mg, Si, O and Fe: the timescale for 
the formation of silicates in the early evolution, when SNae 
enriched the ISM with all of these elements, is therefore 
shorter than the one for carbonaceous grains (carbon is 
also partially locked in the unreactive CO). 

Carbon is injected by SNae (plus a significant con- 
tribution by low-Z AGB stars), but in any case carbon 
is less than the sum of all the other refractory elements 
involved in silicates. Their total number density in the 
ISM, compared to the one of carbon atoms, favors the 
formation of silicates during the early evolution. This can 
be seen in the upper right panel of Fig. \b\ the fraction of 
carbon locked up in CO also plays an important role. In the 
panels of Fig. [7] for the carbonaceous grains and silicates, 
one can note that, as a consequence of the behaviour 
described above, the fraction of silicates in the innermost 
regions grows until the age of about 1 Gyr is reached, 
while the fraction of C-rich dust decreases. Afterwards, 
as AGB stars join SNae to increase the number density of 
Carbon atoms in the ISM, the accretion of carbon in cold 



MCs becomes also efficient, and the silicates/carbonaceous 
grains ratio slowly decreases again. A different behaviour 
takes place in the outermost regions: at the beginning 
there is a ratio similar to the inner regions, but the fraction 
silicates/carbonaceous grains keeps decreasing because (i) 
the accretion in the ISM has long timescales and starts late 
due to the low densities and Z; (ii) the strong contribution 
of C-based dust, due to the low-Z AGB stars, dominates 
in a significant window between SNae and start of the ISM 
accretion. The final result is the high fraction of carbon 
dust with respect to silicates in the outer region. It is 
interesting now to compare our results with the evolution 
of the gradient in Carbon predicted by IZhukovska fe Gail 
(|2009t) . Since in their case SNae are very poor factories 
of silicat es, but significa nt injectors of Carbon-rich dust, 
Zhu kovska fe Gail ( 2009T ) adopt for the early evolutionary 
stages in the inner regions, a mixture almost totally made 
by Carbon-based dust. In our case, on the contrary, we 
obtain in the inner Disk that the partition between the 
two main dust types is more balanced and both play a 
comparable role. In any case, in the outer regions we find 
results very close to those bv IZhukovska fe Gaill (|2009l) . In 
brief, at the current age we all find that the final gradient 
in relative mass fraction of carbonaceous grains has a 
similar positive slop reaching in the outermost region 
a value as high as 70% of the dust mass embedded in 
carbonaceous grains. 

The behaviour of iron dust has the opposite trend: 
it starts very flat and at the current time we have a 
negative gradient, with iron dust significantly contributing 
to the dust mass in the center of the MW. In the early 
stages, we have some contribution to the iron budget by 
SNae (however less than for silicates and carbonaceous 
grains, see top panels in Fig. [5]) and because of the low 
metallicities a negligible contribution by low-Z AGB stars. 
As time goes on, the sum of the contributions from (i) 
CCSNae, (ii) accretion in the ISM that for iron dust starts 
to be significant at t> 1 — 2 Gyr, (iii) AGB stars that reach 
solar metallicities in the inner regions and s tart to inject 
iron d ust, and, finally, (iv) type la SNae (| Piovan et al.l 
l2011al a). The sum of all of these contributions yields a 
gradient with negative slope. 

Radial gradients in element depletion. Once examined 
the results for the radial gradients in element abundances 
both in the gas and dust, we discuss the corresponding 
gradient in element depletion due to dust formation. First, 
we concentrate on the innermost and outermost regions of 
the Disk. In Fig. [8] we show the results for the depletion 
of r C, N, O, Mg, Si, S, Ca and Fe in the ISM in an 
inner ring of the MW Disk centered at 2.3 Kpc. This 
is a region of high density and high star formation in 
turn. The results are shown for both the original (upper 
panel) and c orrected Mg y ields (l ower panel). Indeed, as 
described in IPiovan et al.l (|2011bl ). a small correction to 
the original yields for the Mg under- abundance leads to 
different (better) results for the dust mixture as the Mg 
to Si ratio affects the formation of silicates and the final 
re sults for the SoNe. Th e same group of models presented 
in IPiovan et al.l (|2011bf ) and used to study the depletion 
in the SoNe are shown for the inner and outer regions, 
but limited to those with the coefficient of star formation 
efficiency v = 0.3 and v = 0.7. The infall time scale is 
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varied between 1 and 9 Gyr (or between 3 and 9 Gyr in 
the Mg-corrected models) and three IMFs are examined, 
namely the Kroupa multi power-law IMF, the Larson IMF 
adapted to the SoNe and, finally, the Salpeter one for the 
sake of comparison. We also show the un-depleted (no dust 
formation) abundances of the gas (i.e. the gas+dust total 
abundance). The abundances referring to the same models 
are connected by a line (dashed or continuous for v = 0.3 
and v = 0.7 respectively) to better show how the elemental 
abundances fall down because of the dust formation. As 
expected, all the elements exhibit a lower abundance in the 
gas. Some common features are: (i) models with corrected 
yields of Mg (bottom panels) show a stronger depletion 
because they form more silicates, due to the more balanced 
ratio between Mg and Si; (ii) Mg is much more depleted 
than Si because of the higher abundance A[X gas /H] of Si 
(actually A[X gas /_ff] of Si and Mg are nearly the same but 
the total starting abundance of Si is higher) ; ( iii) compared 
to the depletion in the SoNe examined by iPiovan et al.l 
(|2011bf) . there is a less significant difference between the 
models with v = 0.3 and v = 0.7, with the exception of the 
case with r = 1 Gyr. The reason for this can be tracked 
in the fast onset of the ISM accretion phase in both cases. 
Very early on, accretion drives the dust budget and tends 
to smooth differences due to injections of Stardust; (iv) 
the IMFs with higher number of SNae produce a stronger 
enrichment than the case with a smaller population of 
SNae, but the A[X gas /H] caused by the depletion is quite 
similar. 

The situation for a companion outermost region of low 
density and low star formation centered at 15 kpc is pre- 
sented in Fig. IH1 The meaning of the symbols is the same 
as before. In this Figure we applied a small to the right 
for the case v = 0.7 to avoid superposition of the models. 
Comparing the data displayed in Fig. [5] and Fig. [§] we note 
that: (i) in general the innermost regions form more dust 
with a higher reduction of the gas abundance A [X gas /_£/]. 
This is the result of forming dust in a medium richer in met- 
als and denser in which the accretion process is favored; (ii) 
the effect of the different efficiencies of star formation (pa- 
rameterized by v) is more evident in the low density /low 
SFR regions because it makes more evident the effect of 
metals and ISM enrichment on dust formation; (iii) in the 
outermost regions almost all the elements suffer a compa- 
rable depletion A[X gas /H], simply because the accretion 
process is very weak over there. 

Finally, in the six panels of Fig. [TU]we show the abun- 
dances of [Mg/H], [Fe/H], [C/H], [Si/H], [S/H] and [Ca/H] 
in the gas as a function of the galacto-centric distance. To 
single out the difference between models with and without 
dust depletion, in all the diagrams we show the radial gra- 
dient at the current age for the ISM as a whole (that is 
gas+dust, left sub-panels in each plot) and for the gaseous 
component alone (right sub-panels in each plot) including 
the effect of depletion. Twelve combinations of the parame- 
ters IMF, v and r are shown for both the un-depleted (left 
panels) and depleted (right panels) cases. For each group 
of four panels we have therefore six un-depleted cases on 
the top left sub-panel and six on the bottom left sub-panel. 
The same six plus six cases are shown for the depleted gas 
in the top and bottom right sub-panels. As expected, dust 
globally lowers the amount of elements in the gas across 
the whole disk, but also in some cases alters significantly 



what would be the slope of the gradient in absence of dust. 
In particular, in the innermost regions, where a higher den- 
sity of heavy atoms is available to accrete on the seeds, 
the regularity of the gradient breaks more strikingly than 
in the outermost regions where the gradient keeps roughly 
the original slope. The innermost regions, that are also af- 
fected by the bar (r < 5 kpc) , present a even more irregular 
behaviour, in particular for the elements most depleted and 
participating to the complex process of silicates formation, 
like Magnesium and Silicon. Indeed, Carbon keeps a regu- 
lar gradient as in our model no other heavy elements are 
supposed to participate to the accretion process of carbona- 
ceous grains. Iron is also very irregular in the innermost re- 
gions, since it participates both to the formation of silicates 
and iron dust grains. 

6. Discussion and conclusions 



This i s the last paper of a series of three (|Piovan et al.l 
l2011al fbl) devoted to the study of star-dust and dust accre- 
tion in the ISM. We have analyzed in some detail the chem- 
ical gradients (gas and dust) across the Disk of the MW. 
The MW as a whole and its Disk and SoNe have been con- 
sidered as ideal laboratories in which theories of chemical 
evolution and dust formation are compared with observa- 
tional data. A suitable description of the Disk including 
radial flows of matter and the presence of a inner Bar has 
been adopted to simultaneously explain the bump on the 
inner gas density and the gradients in chemical abundances. 
This model provides the frame to the formation and evo- 
lution of the dust, one of the components of the ISM (the 
gas elemental species being the others). 

The main conclusions of this study can be summarized 
as follows: 

— For pla usible combinations of the parameters, the same 
that in IPiovan et al.l (|2011bh allowed to nicely repro- 
duce the properties of the SoNe, both for gas and dust, 
we find a good general agreement with the mass gradi- 
ents observed in the MW Disk, both for the total gas 
mass and the single element abundances. These latter 
have been compared with the abundances observed in 
different types of objects, like Cepheids, open clusters, 
HII regions, O and B stars and Red Giants. Some strik- 
ing differences between models and observations for C 
and N can be easily explained by means of the stellar 
evolution theories. 

— The evolutionary scheme for dust formation and evolu- 
tion across the Disk, emerging from the simulations, is 
as follows: SNae dominate the dust budget in the early 
stages in the whole disk. As time goes on, AGB stars 
enter the game: they contribute to the C-based dust 
budget if the metallicity keeps low, switch to O-rich 
dust as soon as the metallicity is near solar or super- 
solar. For this reason AGB stars behave in a different 
way in the inner regions where the metallicity quickly 
increases and in the outer regions where the metallicity 
keeps low for long time. The ISM accretion is the main 
contributor to the dust budget across the whole disk. 
It starts very early on in the innermost regions due to 
the faster enrichment and the highest number densities 
of the elements involved in dust and much later in the 
outer regions where metallicity and density are lower 
and where AGB stars of low-Z can be very important 
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for many Gyr. It is clear from this trend that a cru- 
cial role is played by the SNae: the faster they enrich 
the ISM in metals, the faster the ISM starts to play 
a decisive role. Furthermore, the higher the condensa- 
tion efficiencies, the higher is the amount of dust that 
can be injected during the time interval before the ISM 
dominates. The uncertainties on the efficiency of dust 
condensation in SNae ([Piovan et al.ll2011al ) could clearly 
change the pictur e duri ng the early evolutionary phases 
(|Zhukovska et al.ll2f)08l: ICall et al.ll2011allbh . 

— The radial distribution of the dust mass follows the be- 
haviour of the gas showing a bump in the inner zone and 
a negative slope at the present time. It decreases mov- 
ing outwards. The fractional radial mass of silicates has 
also a negative slope, whereas the carbonaceous grains 
have a positive slope. They also play a more and more 
important role as we move toward the outskirts. This is 
ultimately due to the combined effect of several causes 
: (i) the yield from SNae that enrich the ISM with ele- 
ments involved in silicates and dust where star forma- 
tion is more active; (ii) the long time interval during 
which low-Z AGB stars inject carbon in the outskirts 
of the Disk, and (iii) the delay in the onset of the ISM 
accretion in the outer regions. 

— The depletion of the elements is more significant in the 
inner regions, with the highest star formation and den- 
sities, while in the outskirts we have less depletion due 
to lower efficiency of dust accretion and star-dust injec- 
tion. 

— As expected, the depletion significantly alters the radial 
gradients in the gas abundances, in particular for the 
elements involved in silicates, for which (at least in our 
simplified model) many elements enter the process and 
the key element driving the process can vary with time, 
making the whole thing more irregular. 



We can conclude by saying that the picture drawn by 
the chemical model for the Galactic Disk in presence of 
dust is consistent with the way in which SNae, AGB stars 
of different Z and ISM contribute to the total dust budget. 
In particular, since the radial properties of the elements 
along the disk are satisfactorily reproduced, we are confi- 
dent that the dust enrichment scenario described is based 
upon a consistent general enrichment in metals of the ISM. 
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available ones for that element, in such a way to examine the gradients from different sources. The following models are 
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Fig. 5. Temporal evolution of the radial contribution to the abundance of dust by the four types of grain among which 
we distributed the single elements (tha t is s ilicates, carbonaceou s grains, iron dust grains and generic grains embedding 
Ca, N or S - see lZhukovska et all (|2008f ) and lPiovan et al.l (|2011bl) for more details) and the three sources of dust, namely 
SNae, AGB stars and accretion in the ISM. All the contributions have been properly corrected for the destruction of dust. 
Five ages are represented from the early stages to the current time. Upper left panel: the radial contribution of the 
four types of dust grains to the total dust budget at 0.1 Gyr. We show: silicates (continuous lines), carbonaceous grains 
(dashed lines), iron dust (dotted lines) and, finally, generic grains containing S, Ca and N (dot-dashed lines). For each 
group we distinguish the net contributions from the ISM accretion, AGB and SNae, that is: ISM-C, AGB-C and SNae-C 
for carbonaceous grains, ISM-Sil, AGB-Sil and SNae-Sil for silicates, ISM-Fe, AGB-Fe and SNae-Fe for the iron dust and 
finally, ISM-Ca/S/N, AGB-Ca/S/N and SNae-Ca/S/N for the other grains. In all cases the abundance of CO is fixed to 
Ceo = 0.30 taken as the reference value. For the last two ages we also show for the ISM the results for £co = 0.15 and 
Ceo = 0.45. Upper Central panel: the same as in the upper left panel but for 0.5 Gyr. Upper Right panel: the 
same as in the upper left panel but for 1.0 Gyr. Lower Left panel: the same as in the upper left panel but for 5.0 Gyr. 
Lower Right panel: the same as in the upper left panel but for 12.8 Gyr. Since the ratio (Jd/vh grows at varying age, 
the scale of the y-axis is continuously shifted according to the value of the plotted data. 
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Fig. 6. Temporal evolution of the logarithmic radial abundance of dust of (rk,t), in [Mqpc~ 2 ], for all the elements 
belonging to our set and involved into the process of dust formation, that is C, N, O, Mg, Si, S, Ca and Fe (see also 
I Piovan et al.l (|2011bj ) for more details) . All the contributions have been properly corrected for the dust destruction. Nine 
ages are represented from the early stages to the current time, that is 0.1, 0.2 , 0.5, 1, 2.5, 5, 7. 10 and 12.8 Gyr, assuming 
that the formation of the MW started when the Universe was ^0.9 Gyr old (jGail et al.l l2009). Upper panels: from left 
to right the time evolution of the radial abundance o D for Mg and Fe. Central panels: the same as in the upper panels 
but for C, N and Si, from left to right. Lower panels: the same as in the central panels, but for S, Ca and O. Since the 
abundance of changes at varying the age in a different range for each element, the scale of the y-axis is continuously 
shifted according to the represented range of values. 
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Fig. 7. Temporal evolution of the radial mass fraction of single elements or grain families embedded into dust 
erf (rfc , t) /cr D (r*fc , t) , normalized to the total dust budget, for some of the elements belonging to our set and involved 
in the process of dust formation (that i s C, N, O, Mg Si, S, Ca and Fe) and some grain types (silicates, carbonaceous 
grains, iron grains and Ca/N/S - see also lPiovan et al.l()2011b[ ) for more details). All the contributions have been properly 
corrected for the dust destruction. Nine ages are represented as in Fig. El from the early stages to the current time, that 
is 0.1, 0. 2, 0.5, 1, 2.5, 5, 7 , 10 and 12.8 Gyr, assuming that the formation of the MW started when the Universe was ~0.9 
Gyr old fG ail et al.ll2009f) . Upper panels: from left to right the time evolution of the radial mass fraction of Mg and 
Fe, normalized to the total dust mass. Central panels: the same as in the upper panels but for O and Si, from left to 
right. Lower panels: the same as in the central panels, but for S, N and Ca (lumped together) and O, silicates (that is 
the mass fraction of Mg, Si, Fe and O involved into quartz/pyroxenes/olivines) and C (that in practice, once subtracted 
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Fig. 9. Depletion of the elements C, N, O, Mg, Si, S, Ca and Fe in the ISM for an outer ring of the MW at 15.0 Kpc. 
Top panels: the element depletions are obtained using the original yields, no correction for the Mg under- abundance is 
applied. The choice of the parameters and the meaning of the symbols are the same as for Fig. [5J For the sake of a better 
representation of the data, the case with v = 0.7 is slightly shifted to the right in such a way to avoid the superposition 
with v = 0.3. As in Fig. [51 the dashed lines connect the un-depleted and depleted gas abundances for the v — 0.3 case. 
Bottom panels: the same of the top panel, but with the original yield corrected for the Mg under- abundance. The same 
timescales of Fig. [8] are used. 
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Fig. 10. The Galactic abundances of [Mg/H], [Fe/H], [C/H], [Si/H], [S/H] and [Ca/H] in the gas as a function of the 
galacto-centric distance. In all diagrams we present the radial gradient at the current age for both the ISM component 
as a whole (that is gas+dust, left panels in each plot) and the gaseous component alone (right panels in each plot) 
taking into account the depletion of the elements due to the presence of dust. Twelve combinations of the parameters 
are considered for both the un-depleted (left panels) and depleted (right panels) cases. The legend has been split between 
left and right panels for the sake of clarity. The following cases are shown. Kroupa IMF: (t, v) — (3, 0.3) (continuous 
green line), (t,u) — (6,0.7) (continuous blue line), (t,v) = (9,0.3) (continuous red line) and (r, v) — (9,0.7) (continuous 
blue line); Salpeter IMF: (t,u) = (3,0.3) (dotted green line), (t,v) = (6,0.3) (dotted blue line), (t,v) = (9,0.3) (dotted 
red line), and (t,v) — (9,0.3) (dotted blue line); Larson SoNe IMF: (r, v) — (3,0.3) (dashed green line), {t,v) = (6,0.3) 
(dashed blue line ), (t,v) = (9,0.7) Gy r (dashed green line) and (r, v) — (8,0.3) (dashed blue line). All theses cases are 
in lPiovanet"aTI(|2011bf) for the SoNe. 
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